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i Abstract 
a = MS 
q The complex multifaceted decision situations present today suggest 
iq the need for a timely, automated too. A-.decision-maker is for ed into 
3 comparing alternative actions or systems over an entire set of different 
a measures of merit. This effort is an on-line, real-time, comp ter-based 
a decision aid designed to assist the decision-maker in clarifying prefer- 
4 ences in a complex decision environment. It is applicable to problems 
4 which may be represented by a hierarchy of objectives to be satisfied. 
4 The program is MADAM: Multiple-Attribute Decision Analysis Model, and it 
a is written in FORTRAN V and is implemented on the CYBER 175 system. MADAM 
4 is designed to aid the decision-maker as he or she progresses through prob- 
a lem formulation, parameterization, sensitvity analyses, and a decision, 
: including storage of all data and rationales. Deterministic problems are 
3 analyzed through Multi-Attribute Utility Theory concepts and an additive 
3 value function is utilized for sensitivity analysis. Pairwise preferential 
4 independence is tested between AL EnVubeS sf The sensitivity analyses in- 
P clude a cumulative weight anaiysis, a relative weight analysis, and an 
¥ attribute level analysis. The analyses may be conducted by fixing an ob- 
a jective to be considered and conducting the analysis across the alternative 


| 
{ systems or actions, or conversely by fixing the alternative to be consid- 
2 | ered and conducting the analysis across a desired set of objectives. 
3 | fine werk is divided into two volumes. Volume I is a theoretical pre- 
{ sentation and includes a user's manual. It requires no programming ex- 
{ pertise anc may be used independently of Yolume II, Volume If is a pro- 


gramming manual including the source code. 


\" nay not be used independently 


| 


par 
Re | i of Volume I. 
| 


I INTRODUCTION 


This manual is designed to supplement the information contained in 
Volume I of this thesis. The foundamental intention to this volume is to 
provide decision analysts with an elucidation of the details of MADAM in 
order to better understand the program and, hence, facilitate its imple- 
mentation. Future analysts may also have the desire to modify or extend 
the scope of problems suitable for analysis with MADAM, and the following 
information will be very useful. 

MADAM is written in FORTRAN V and implemented on the CYBER 175 system, 
An effort was made to keep the code as machine independent as possible, 
and for those situations where this was not possible, an explicit discus- 
sion concerning modificetions pre-requisite to transporting MADAM IS pro- 
vided. All references to an example problem will refer to the example 


iVustrated in Chapter V of Volume I. 
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II DATA STRUCTURES 


MADAM utilizes a random access file for storing and manipulating ‘data. 
Throughout the program are pairs of, subroutines and functions which are re- 
sponsible for setting and retrieving data, This is possible because FORTRAN 
V will not recognize an array variable until it is determined that ho func- 
tion by that name exists. Thus, for example, MADAM uses the function ARAY 
(I, J) to represent a pseudo-array variable. A complete list of the sub- 
routine function pairs which represent psuedo-array yartables is given in 
Table 1. Each record of the random access file is 100 words in length. 
This 100 word record is divided up in to data fields as shown in Figure j 
for the first 500 records. The remaining records have data fields as shown 
in Table 2. 

NOTE***The current ability to access both. real and 
integer (character) data from a record is dependent 
on the fact that the CYBER 175 system uses one word 
for both types of information. In order to imple- 
ment MADAM on a system which uses two wards for real 
data and one word fur integer data, it would be nec~ 
essary te modify the psuedo-array yartables so that 
the desired words will be accessed. This will re- 
quire manipulation of the record length and storage 
(retrieval) subscripts. 

Each data file (tree) constructed will require about 300 PRU for per- 
manent storage. 

As illustrated in Figure 1, each node record contains information 


about the position of the node in the objective hierarchy plus any node- 
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Table 1 Pseudo-array Variables 
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Node Referance Number Digtt 


Data 


First Descendant 


First Crosslink 


Backlink 


Relative Weight 


Cumulative Weight 
Associated Attribute 
System Values 
Comments (Rationale) 


| Objective 
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Figure 1 Data Fields for Node Recard 
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Node Record 


Word(s) 


8 to (NSYS+7) 
(NSYS+7) te 64 


65 to 80 


Number of Nodes (NNODES) 
Number of Attributes (NATT) 
Number of Systems (NSYS) 


Null Record 


Regression Intercepts 


Regression Slopes 


Sum of Squared Errors 


Regression Form Indicators 


System Names 
Attribute Names 


Null Record 


Worst Attribute Levels 
Best Attribute Levels 


wena 


fable 2 Record Content 


Variable(s) 


-none- 
OBUCTV (1) 


to OBUCTY (16 


Record No. 


1 to 500 
501 
502 
503 
504 
505 
906 


ae 
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\ specific data. The positional information is provided by the first four 
words in the record. These words contain the node reference number digit, 


the first descendant, the first crosslink, and the backlink respectively. 
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The iode reference number digit indicates where this node is on its span. 


For instance, in the sample problem, the node reference number digit of 
the node with the objective: to maximize aircraft survivability, is two. 
The first descendant pointer contains a number of the first child of the 
given node (if any). Absence of a first child (IRAY (1,2) < 0) impltes 


that the given node is a data node. The first crosslink pointer contains 


a number of the immediate siblings of the given node (if any). The back- 
link pointer contains the number of the node which has this given node as 
a first descendant or a first crosslink. 

The remaining words in a node record provide the data for identifying 
the node and incorporating it into the decision analysis. The fifth word 
of the record (ARRAY (1,1)) contains the weight of the given node relative 
to its siblings. These relative weights are rormalized to a unity sum 
across the span. The sixth word of a node record (ARRAY (1,2)) contains 
the cumulative weight of a given node. It is automatically calculated by 
MADAM during the "folding-back" of the objective hierarchy (see subroutine 
CALC), and it is a function of the node's relative weight and the cumulative 
weight of its parent. The seventh word in a node record (LABEL (1)) con- 
tains the attribute label to be associated with that node, Jt is relevant 
only for data nodes, and is eft a null field for all nodes interior to the 
tree. 

NOTE***A1) character data is designed around the 
CYBER's 10 letters/word ferinat. In order to 


implement MADAM on s system with a different 


storage format, it will be necessary to scale 
all input data or data structures to incorporate 


a different number of characters per word. 


The next set of words in a node record (YVRAY (1,J)) hold the values 


of the alternative systems relative to a given node. These values are 


either calculated by MADAM (for an interior node) or indirectly input by 
the user (for a data node). Only as many words as needed are used for 
this information basea on the number of systems. Those words between 

the system vaiues and the 64th word are used to hold comments or ration- 
ale concerning either the relative weighting of the children of a given 
Thus, the 


node or the input attribute levels of the alternative systems. 


amount of space for these comments is inversely proportional to the number 


. 

a 

| of alternative systems. Also, because of this structure, the number of 
F 

I ( alternative systems is limited to 57, Finally, the objective associated 
with the given nede is contained in the 64th to the 80th words. 
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a c Ul SENSITIVITY ANALYSIS 


~ The sensitivity analyses are controlled by tha routines in the pri- 
mary overlay (4,0). The types of analysis handled by the program are: 
relative weight analysis (RSA), cumulative weight analysis (CSA), attri- 


bute level analysis (LSA), and a system analysis (SSA) based on one of 


the first three tvpes. The first three types of analysis involve deter- 
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mining a fixed node and conducting the appropriate analysis for all sys- 


> i tems at that node. The last type of analysis (actually three subtypes 
io SSA/RSA, SSA/CSA, SSA/LSA) involves fixing the alternative system and then 
- i 
a ; conducting the appropriate analysis over a set of nodes, 
{ 


The entry pcint of the overlay, routine SENSTY is used to determine 
the main type of analysis (CSA, RSA, LSA, or SSA). For any of the first 
three types of sensitivity analysis, control is shifted to routine DETNOD 
which is used to fix the analysis to the desired node, Routine DETNOD also 


establishes the range over which the factor (cumulative weight, relative 
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weight, or attribute level) is to be considered. For the last type of 
analysis (SSA), control shifts to routine SYSSEN which in turn calls upon 
the routines, SENVAL, SENTYP, and NCDSET to establish the fixed system, 
the factor to be considered, and the set of nodes over which the analysis 
is to be conducted. 

Regardless of the type of analysis to be performed, once the initial- 


ization process {as described in the preceding paragraph), has been com- 
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pleted, control shifts to the routine COMPUT. This routine utilizes the 
equations presented in Chapter IV (of Volume I) to generate a matrix of 


overall values. Each row of this matrix corresponds to a different level 


‘ { 
N of the factor under consideration. For the CSA, RSA, and LSA, each column 


C corresponds to a different system. For the SSA, each column represents a 


different node. : a 


After the overall values haye been calculated, the routines GRAPH 
and TABDIS are used to display the information in a graphical or tabular. 
manner respectively (see Figures 4.1 through 4.4 in Volume I}, At the 
conclusions of the desired type of sensitivity analysis, control jis eupned 


over to the main option selection. 
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TV HIERARCHY MANIPULATION 
Subroutines NODIN and FINO 

These subroutines are used to allow the user to enter the tree struc- 
cure at a desired node. If the desired node is not an element of the tree, 
FIND will identify that node which is closest to the input node. All nodes 
are referenced by their node reference number (NRN). Except for option 
***MOD***, the input of a non-validated NRN will be interpreted as an in- 
dication tnat the user desires to exit from the option, Program contro] 
will be return to the calling portion of MADAM, For option ***MOD*** an 
invalide NRN is interpreted as a desire to introduce a new node having 
the input NRN, so the appropriate branches (and dummy connecting nodes) 
are constructed so as to create a valid environment for the input node, 

Subroutine NODIN. This routine allows the user to input a desired 
NRN. If anything but a null string (space followed by <CR>) is entered, 
this routine will call on FINO to establish the existence of a node cor- 
responding to the input NRN. 

Subroutine SIND. This routine implements a modified breadth-first 
search of the tree to search for a node corresponding to the input NRN, 
With each digit of the NRN representing a level, FIND starts with the 
first digit and attempts to locate a node with the desired NRN digit by 
crossing the span at that level. This process is repeated for all re- 
maining levels (digits) of the input NRN. If a match cannot be estab- 
lished at a given level, this indicates that the desired node would be 
an added crosslink at that level. If a match is established, the next 
digit is brought up for consideration. If no more downlink path is 


available (a dat node has been reached) the routine is terminated. 
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“Figure 2 iliustrates an example of this search process. 
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No Match 


No Match 


Second digit is satisfied 


No Match 


All digits satisfied so 
exit 


Figure 2 The Search Process (FIND) 
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These routines driye a depth-first search of the tree structure, Sub 
routines PRENEX and PRETAT are used to initialize the searcli process at a 
desired node or at the root node respectively. Once all variables have 
been initialized, subroutine NEXT performs a depth-first traversal, advanc- 
ing one node each time it is called. J 

Subroutine PRENEX. This routine elicits a NRN from the user by ine 
voking NODIN. If a valid NRN is entered, the traversal wilt begin at the 
node corresponding to the input NRN. All descendants of this input node . 
will be encountered in the resulting transversal. 

Subroutine PRETOT. This routine initialized the variables required 
by NEXT so as to run a transversal of the entire tree structure, It is 
usually invoked just prior to a potential application of PRENEX, thus al< 
lowing the user to gracefully opt a tour of the entire hierarchical struc- 
ture. 

Subroutine NEXT. This routine tours the hierarchical structure begin- 
ning at the node specified by either PRETOT or PRENEX. Since a check is 
made as to the validity of the input node before any processing is done, 
it is pessible to exit the current option by entering a null string during 


option PRENEX. The logic implemented in NEXT is illustrated in Figure 3. 
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Figure 3. Logic Flow for NEXT 
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Subroutiné SPAN 


This routine sions for the input of the objective hierarchy struc- 
ture by sequentially providing nodes in a depth-first fashion. The user 
inputs immediate descendants of the current node, Subroutine CHECK is _ 
used to establish a satisfactory set of children for the current node. 3 
‘Once a satisfactory set is established the first descendant of the 


current node becomes the new parent node. When a data node is reached 


PE eee nen ee Ty 


(no descendants are input); the crosslink to the current node becomes 
the new parent node. If no more siblings exist, the routine wil) back 
up one level and establish the first crosslink as the new ;parent node 
(see Figure in the main thesis text). 

This routine calculates intermediate values (collapses the. tree) 
by implementing a modified depth-first search. The tree is examined by 
each path from the root node to a data rode, and the cumulative weights 
are Calculated. System values based on the attributes are then computed 
for each node. The logic flow of CALC is provided in Figure 4, The | 
cumulative weight of each node is given by: 

CUMNT (NODE) = CUMHT(PARENT)* RELWT (NODE) 
The system value at a particular node are the products of the system 


value at a data node and the cumulative weight of the upper nodes for 


each system. 
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Figure 4 Logic Flow of CALC 
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Variables Used 


ne ee : a 


The followitig i's a list of variables used in the program. This list 


includes only important variables and does not include scratch variables, 


do-loop indices, arid dummy arguments. Psuedo-array variables are in the 


next section. 


AIN 
ARRAY 


ATTMAX 


ATTMIN 
ATTNAM 


N 


AT?RI 
ATTR2 
BO 

BOEXP 


BOLN 
BOROOT 


BOSQR 
BY 


BIEXP - 
BILN 
BIROOT - 
BISGR - 
BAND - 


CMD ~ 


vector of relative weights 


used to vead and write to the data file. It holds one 
record worth of information 


maximum level or attribute 
minimum level of attribute 


name of attrihute being perturbed during sensitivity 
analysis 


name of first attribute 

name of second attribute 

intercept of linear regression (linear form) 
intercept of linear regression (exponential form) 


intercept of linear regression (logarithmic form) 


~ intercept of linear regression (square-root forin) 


intercept of linear regression (squared form) 
slope coefficient of linear regression (linear form) 
slepe coefficient of linear regression (exponential form) 


slope coefficient of linear regression 


( 
( 

slope coefficient of linear regression (logarithmic form) 
(square-root form) 
( 


slope coefficient of linear regression (squared form) 
boundaries of indifference curve window 


main option selection 
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DELI 
DEL2 
ERRELG 


FACLVL 
FINDST 
ICHECK 
TCONT 
IDATA 
TFADD 
IFIND 
INDERR 
INDX 
INDXST 
IQUIT 
IRRAY 


ISAYD 


ISIN 


~ 


change in first attribute :: 
change in second attribute 


counts the number of consecutive errors in user entry 
during main option selection 


level of attribute 
set of IFIND for SSA 


flag if using between-node check 
0 - not using check 
1 - using check 


flags whether or not to continue processing 
0 - stop 
1 - continue 


flags whether or not node is data node 
Q - interior node 
1 - data node 


flags relationship of current node to input node. 
Determines how to add th input branch to the existing 
branch. 

2 - IFIND is parent to input node 

3 - IFIND is sibling to input node 


points to tne last node of the branch contained in INNRN 
flaq for PPI 

0 - PPI hold 

1} - PPI fails 

NRN digit for IFIND 

set of INDX for SSA 

Tlag to continue processing 

0 - continue 

1 - stop 


used to read and write to the data file. It holds one 
record warth of information. 


flags whether current information is saved on data file 
in case of abort 

0 - not saved 

1 - saved 


record number of record currently in memory (in ARRAY or 
IRRAY ) 
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ISYS - number of systems being considéred during sensitivity 
analysis . 

ITOL. - tolerance during PPI testing 

ITOTL = - flags type of tree traversal 


0 - not a total tree traversal 
1 - traverse all -descendants 


KIDSET - set of NCHILD for nodes during SSA 

LEVEL - vector containing branch. to node IFIND 
(*,1) - location of record containing node * 
(*,2) - flags whether or not visited by PRETOT-PRENEX-NEXT 
(*,3) - used in copying 

LINE ~ one line of data for plot 

LVL - length of LEVEL vector 

MATRIX - pairwise-comparison matrix 


MAX? maximum level of first attribute 
MAX2 


MINI 


maximum level of second attribute 


minimum level of first attribute 


MIN2 minimum level of second attribute 


MTYPE type of SSA to be done 


NATT nuroer of attributes 


t 


NCHILD number of children for a node 

NCROSS - number of siblings on a span 

NDEEP - depth of tree (maximum number of levels) 

NDIFF - number of digits in INNRN not matched by routine’ IND 
NLOUD - number of allowed nodes plus four 

NLVLS  ~- Tength of INNRN 

NNODES - number in tree plus one (master) 

NSYS - number of alternative systems 


NTAPE - defines which data files is in use 


NTREE  - maximum number of simultaneous data files 


C4 .6RdGtY ~ objective associated with ‘a-particular node 
SENS = type of sensitivity analysis: ~~ 
SSE - sum-of squared error (Tinear form) 
SSEEXP - sum of squared error (exponential form) - = °° 
SSELN  ~ ‘sum:of squared error (logarithmic form) 
SSERT - sum of squared error (square-root form) 
SSESQR - sum of squared error (squared form) 
TOLER - tolerance for PPI testing 
TOP ~ objective of parent node during between-node check 


USER - user identifier 


VAL attribute levels which dafine data points for the 


individual value function 


WHOLD - matrix of overall value of systems during sensitivity 
analysis 


(G WMAX = maximum level of factor during sensitivity analysis 


WMIN minimum level of factor during sensitivity analysis 


Psuedo-arry Variables 
The following are functions which are treated as if they were array 


variables (s3e Table 1). 


ARAY(1,0) relative and cumulative weight of node 
. ATT(T,d) best/worst levels of attributes 
ATTU(T) attribute names 
CRAY(T) comments or rationale (subroutine) 
IRAY(I,J) node position information 
LABEL (I) name of attribute associated with node I 


OBJECT(I ,OBUCTY) objective at node | 


PARAN(T ,J) regression parameters 
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SYSLBL(T) alternative: system names 

VRAY(1,J) system values at node I 
Cross Reference Map 

The listing contains a partial load: map, and combined with the com- 
ments within each program unit, a corss reference map of where each vari- 


able is used or modified may be generated. 


CX VI. Program Structure 


MADAM is divided up into five overlays (one main, four primary) 
which reduce the required field length for loading to about 60K. These 
overlays are machine dependent, but they may be removed by deleting the 
overlay statenents, changing the program statement following each over~ 
lay statement to a subroutine statement (except for PROGRAM MADAM), and 
changing all overiay calls to subroutine calls to the subroutine making 
entry into the described overlay. Each primary overlay is designed to 
operate independently of the other primary overlays. Only the main over- 
lay and one primary overlay are in the user's field length at one time, 
thus considerably reducing the required core memory allocation. The 
following is a list of routines in MADAM in the order in which they ap- 


pear in the source code. 


OVERLAY (XFILE,0,0) - main overlay 
PROGRAM MADAM 
FUNCTION ARAY(T,d) 
SUBROUTINE ASET(I,J,V) 
FUNCTION ATT(I,3) 
FUNCTION ATT1(1) 
SUBROUTINE ATTLBL(1 ATTRIB) 
SUBROUTINE ATTSET(2 ATTRIB) 
SUBROUTINE ATTVAL(I,J,V)} 
SUBROUTINE COPYR 
SUBROUTINE CRAY(T) 
SUBROUTINE CSET(I) 
SUBROUTINE FIND 
SUBROUTINE HELP 
SUBROUTINE INTRO 
FUNCTION IRAY(I,J) 
SUBROUTINE ISET(I,S,M) 
FUNCTION LABEL(I) 
SUBROUTINE LSET(I,LABELT) 
SUBROUTINE NEXT 
SUBROUTINE NODIN 
SUBROUTINE ORJECT(I ,OBUCTY) 
SUBROUTINE OSET(I ,OBUCTY) 
FUNCTION PRRAM(I,J) 
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SUBROUTINE PRENEX 
SUBROUTINE PRETOT 
SUBROUTINE PSET(I,J,V) 
FUNCTION SYSLBL(1) 
SUBROUTINE SYSSET(I,SYSLBL) 
SUBROUTINE TLOAD 
SUBROUTINE TSAVE 
FUNCTION VALU(X LABEL) 
FUNCTION VRAY(I,0) 
SUBROUTINE VSET(I,K.V) 
FUNCTION XLEVEL(X, LABEL) 


OVERLAY (XFILE,1,0) - primary overlay 
PROGRAM DUMMY 

SUBROUTINE ADDSYS 
SUBROUTINE ASK(I, XV ,X2,XMID) 
SUBROUTINE CHECK 
SUBROUTINE DELSYS 
SUBROUTINE GRAFIX 
SUBROUTINE INITT 
SUBROUTINE MASTER 
SUBROUTINE MODIFY 
SUBROUTINE NEW 
SUBROUTINE PICTUR( INDEX } 
SUBROUTINE PPI 

s SUBROUTINE PRUNE 
SUBROUTINE ROATT 
SUBROUTINE RGSYSL 
SUBROUTINE RDTTL 
SUBROUTINE REGRS( INDEX) 
SUBROUTINE SPAN 
SUBROUTINE STAT 
SUBROUTINE TAPER 
SUBROUTINE VALUE 


OVERLAY (XFILE,2,0) - primary overlay 
PROGRAM WYLOAD 

SUBROUTINE CALL 

SUBROUTINE DVIDE(N,A,D) 

SUBROUTINE CUVEC(NCROSS , [FIND AIN} 

SUBROUTINE ROV 

SUBROUTINE RDWT 

FUNCTION SuM(N A) 

SUBROUTINE NYLODI (OPT) 
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QVERLAY (XFILE,3,9) - primary overlay 
PROGRAM NUM 

SUBROUTINE DISPLA 

SUBROUTINE DSPLOT 

SUBROUTINE NEWPG 

SUBROUTINE NUMREV 


we oe se & 


o OVERLAY (XFILE,4,0) - primary overlay 
PROGRAM SENSTV 
2 SUBROUTINE COMPUT 
SUBROUTINE DETNOD 
SUBROUTINE HEADER(I ,OBUCTV , ANSWER} 
SUBROUTINE SENVAL 
SUBROUTINE TABDIS 
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